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Abstract 

Non-uniform black strings coupled to a gauge field are constructed by a perturbative method 
in a wide range of spacetime dimensions. At the linear order of perturbations, we see that the 
Gregory-Laflamme instability vanishes at the point where the background solution becomes ther- 
modynamically stable. The emergence/vanishing of the static mode resembles phase transitions, 
and in fact we find that its critical exponent is nearly 1/2, which means a second-order transition. 
By employing higher-order perturbations, the physical properties of the non-uniform black strings 
are investigated in detail. For fixed spacetime dimensions, we find the critical charges at which 
the stability of non-uniform states changes. For some range of charge, non-uniform black strings 
are entropically favored over uniform ones. The gauge charge works as a control parameter that 
controls not only the stability of uniform black strings but also the non-uniform states. In addition, 
we find that for a fixed background charge the uniform state is not necessarily the state carrying 
the largest tension. The phase diagram and a comparison with the critical dimension are also 
discussed. 



1 Introduction 



It is known that black objects with translational symmetries such as black branes and black strings 
suffer from the Gregory-Laflamme (GL) instability, breaking the translational symmetries [DO- An 
expectation is that they decay into individual black holes, having larger entropy. This naive expec- 
tation, however, was called into question because of the result by Horowitz and Maeda [3], which 
demonstrated that a black-string horizon cannot pinch off in finite affine time on the horizon. Al- 
though a dynamical simulation was reported [I], the endpoint of the collapse remains an outstanding 
open problem, and several possibilities have been pointed out [SJE]. 

As a solid approach, there are extensive studies of static black objects in compactified spacetimes, 
such as Kaluza-Klein (KK) spacetimes. First, a non-uniform black string (NUBS) is constructed 
perturbatively in D = 5 (five-dimensional) vacuum spacetime [7J, which was generalized to arbitrary 
dimensions and the critical dimension was discovered [81 O 110] : the order of phase transition from 
uniform to non-uniform black strings is second order for D > 14 (D > 13) in a microcanonical 
(canonical) ensemble, while it is first order for D < 13 {D < 12). Recently, such critical dimensions 
were reproduced by an approach based on the Landau theory of phase transition |11[ [T2] . Regarding 
the localized black holes in KK spacetimes, a perturbation method was developed, and small localized 
black holes were constructed [I3j [HI HU [EJ CE7J . After these perturbative analyses, the phase diagram 
of a black-hole and black-string system in fully non-linear regime was clarified in a series of papers |18[ 
fl~9j [20] (see also [21] for improved numerical analysis for black string branch) . The result suggests that 
the branches of static black hole and black string merge at a topologically changing solution. (See 
recent numerical work [2H I22j.) 

Although the GL instability is inevitable for neutral black strings/branes, charges can prevent 
them from suffering from the instability. On this point, there is an interesting conjecture called the 
Correlated Stability Conjecture (CSC) (or called Gubser-Mitra conjecture [SIGH), which asserts that 
the GL instability of black objects with a non-compact translational invariance occurs iff they are 
(locally) thermodynamically unstable. This conjecture has been supported by several examples [251 
[26] [27] . Although we speculate that the charges (e.g., a charge associated with a gauge field, KK 
momentum, angular momentum, and so on) play crucial roles in perturbative and non-linear regimes, 
our understanding of the stability and phase structure of charged black strings/branes is restricted to 
very special cases |28 [llOU 29j. in which the systems can be translated to a vacuum system. The aim of 
this paper is to see the non-linear effects of a charge which stabilizes the uniform states at the linear 
order, as the CSC asserts. For this aim, we perform the higher-order static perturbations of magnetic 
black strings. And then, the physical properties of constructed solutions are investigated. The main 
result is that for fixed spacetime dimensions there exist critical charges at which the order of phase 
transition changes in the sense that charged NUBSs become stable, compared with the uniform ones. 
In other words, it means that the stable inhomogeneous black strings are allowed to exist in each 
dimension. This is in contrast to the vacuum case, in which such stable states emerge only in higher 
dimensions {D > 14). 

The organization of this paper is as follows. In Sec. [2] we set up the problem by introducing an 
action integral, background solutions, and so on. In Sec. [3] the linear perturbation is analyzed to 
some extent. There, an interesting application of Landau theory of phase transition to charged system 
is discussed. In Sec. S] we describe the higher-order perturbations in detail. With the numerical 
results of Sees. [3] and [H we discuss the physical properties of charged NUBSs in Sec. [5j Section [6] 
is devoted to summary and discussion. Appendix [A] includes the general form of explicit Einstein 
equations. We only consider magnetically charged black strings in this paper. The possible other 
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components, corresponding to an electric field, are discussed in Appendix [Bl The calculation of some 
thermodynamic quantities are described in Appendices O and |Dj 



2 Setup 

2.1 Action and background solution 

We consider the following (d + l)-dimensional action (d > 4): 



1 



2(d-2)! ' 



•2 

r d-2 



where G d +i is a (d+ l)-dimensional gravitational constant and 3^-2 is a (d — 2)-form field. We denote 
the total spacetime dimensions by D = d + 1. The equations of motion (EOMs), obtained by varying 
the action, are 

p _ 1 rp M2---Md-2nr j ~ ^ n-2 

/ " /_ 2(d-3)! " J ^-^-2 2(d-l)! 5 ^ ' 

d li (y/ = g9* f ""' ta - a ) =0. (2) 

In addition, the form field satisfies the Bianchi identity, dGFd-2 = 0. 

In this paper, we find non-uniform black string solutions emerging from the GL critical point. 
Given the fact that the GL critical mode is a s-wave one |30j . it will be convenient to begin with the 
following static axisymmetric metric ansatz: 

ds 2 +i = _ e 2A(r,,)^2 + e 2B(r,z) ^(^2 + + ^(r,^^ (3) 

Here, we assume that the horizon of black string extends along the z-direction. For uniform black 
strings, the functions A, B, C depend only on the radial coordinate r. Although the function H(r) can 
be eliminated by a radial gauge transformation, we keep H(r) for later convenience and generality. 
The Ricci tensor and explicit EOMs for this metric ansatz are included in Appendix [A] With this 
metric ansatz, we can find a general solution of magnetic field by solving Bianchi identities: 

3~ = Q m e d _ 2 , (4) 

where Q m is a constant and Ed-2 is the volume element of a unit (d — 2)-sphere. In the rest of this 
paper, we only consider magnetic solutions. Thus, the general solution ([3]) of the form fields makes 
our analysis very simple; we do not need to consider the perturbations of the form field independently. 
For d = 4 and d = 5, the gauge field can carry an electric charge. In particular, the 3-form field 
becomes self-dual in d = 5 and the electric field gives the same results as those for the magnetic case 
(Appendix |B]). 

By the dimensional reduction method, the uniform black string and black brane solutions of the 
action ([T]) can be constructed from the black hole solutions in a d-dimensional dilatonic system [31J . 
The explicit expression of a magnetically charged black string is given by 

ds 2 d+1 = -f + dt 2 + -r^dr 2 + f.dz 2 + r 2 dn 2 d _ 2 , 
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f±(r) = 1 

<F = y/(d - l)(d - 3) (r+r_)( d - 3 )/ 2 s d _ 2 . (5) 
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KK 
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0(e) 




Xlfi 




k ~ O(l) 


0(6 2 ) 


-^0,1 


-^2,0 






0(6 3 ) 




-^3,0, 




fei ~ 0(e 2 ) 



Table 1: Table of unknown functions up to the 3rd order. A mode X np appears at 0(e n+2p ). The zero 
modes Xq jP appear at 0(e 2p ). In the last column, the wavenumber ko of order O(l) and its corrections 
are listed to show at which order they enter into the decoupled ODEs. 



For < r_ < r+, this solution has an event horizon at r = r + and an inner horizon at r = r_ where 
a curvature singularity exists. Introducing new variables by 



A = a + lny%, B = b + \n^JZ, C = c, H = ~ In [/+/ 2 ] , 



a deformed black string can be described by 



ds 



d+l 



-e 2a ^ f+dt 2 + J*,*) f_ 



dr 2 
1+T- 



+ dz< 



+ 



e 2c ^r 2 dn 2 



d-2- 



(6) 



(7) 



Substituting the transformation ([6]) into the general form of Ricci tensor (|55p and source term (I56j) . 
we obtain EOMs for a, b and c. 



2.2 Perturbation scheme 

Following [7j, we can construct non-uniform black strings/branes by static perturbations, as well as to 
specify the GL critical point quite easily. Before the detailed discussion, which will be given in Sees. [3] 
and[U let us introduce the general perturbation scheme here. 

First, we rescale coordinates r, z, and "charge" parameter r_ by the horizon radius: 



r z 
y = — , x= — 



r + 



(8) 



Then, we expand the metric function X(x,y) (X = a,b,c) around the uniform solution, 

oo 

X(x,y) = ^2e n X n (y)cos(nKx), 



n=0 

oo 



X n (y) = £ e 2 PX n , p (y), K = J2 ^ 

p=0 q=0 



(9) 



where Xo t o(y) = is imposed. Here K is the GL critical wavenumber, in other words, L = 2n/K 
gives the asymptotic length of the compactified space, and e is an expansion parameter. Substituting 
these expansions into the Einstein equations, we obtain ordinary differential equations (ODEs), to be 
solved order- by-order (Tabled]). 

To have an insight into the whole perturbation analysis, it would be nice to see the linear per- 
turbations. For simplicity, let us focus on d = 5, and we always restrict ourselves to this number 
of dimensions when we write down perturbation equations hereafter. The structures of EOMs and 
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boundary conditions are similar for the other dimensions (see Eqs. (|55p and (|56p ). At the linear order 
0(e), we have the X\ t Q mode. By combining Einstein equations, we obtain coupled ODEs, which 
determines 01,0(2/) an d 01,0(2/): 

£Si ] oi l0 + P [1] ci,o = 0, 

L®c lfi + pWa 1 ,o = 0, (10) 
where l)£ and P$ (i = 1, 2) are linear operators given by 

f [i] = 62/ 4 - 92/ 6 +g 2 (-8 + 82/ 2 + 32/ 4 ) d 4g 4 - 4g 2 + A; 2 2/ 6 (2 _ 3y 2) 



dy 2 2/(g 2 - 2/ 2 )(2 - 52/ 2 + 32y 4 ) dy (q 2 - 2/ 2 ) 2 (2 - 5y 2 + 32/ 4 ) 



£[2] = _oP_ 2[-42/ 2 + 62/ 4 - 32/ 6 + g 2 (8 - 13y 2 + 62/ 4 )] rf 
fe ° dy 2 y{q 2 -y 2 ){2-5y 2 + 3y*) dy 

12q 4 + Ay 4 + 2At§^ 6 - 3fc§7/ 8 + 4q 2 (3 - Wy 2 + 3y 4 ) 

3[22/ 2 - 32/ 4 + g 2 (-6 + 72/ 2 )] d I2q 2 (l + q 2 - 2y 2 ) 



pm 



2/(g2-2/2 )(2 _ 5y2 + 3 y 4) dy (q 2 - y 2 f{2 - 5y 2 + 3y*) ' 



S[ 2] _ (-1 + ?/ 2 )[3g 2 (-2 + y 2 ) + 2/ 2 (2 + y 2 )] d 4[g 4 + 2/ 4 -g 2 (l + 2y 4 )] nn 

y (g2_ y 2 )(2 _ 5y2 + 3y 4 ) dy {q 2_ y 2 ) 2 {2 _ 5y 2 + 3y 4y I ^ 

The function 61,0(2/) is n °t independent but can be written in terms of 01,0(2/) an d Ci 0(2/), 

, , , (-1 + g 2 )y 2 ai,o + (-1 + 2/ 2 )[(6g 2 - 3y 2 )c lfi + y(q 2 - y 2 )(a lfi + 3ci, )] 

&1 ' 0(2/) = (g 2 -2/ 2 )(-2 + 32/ 2 ) • (12) 

As we will discuss in Sec. [31 Eq. (|10p can be solved as a shooting problem. The equations for higher- 
order perturbations can be reduced to similar systems of ODEs. 

2.3 Physical quantities 

Here, we prepare the expressions for physical quantities of the metric ([7]), with which we will later 
discuss the thermodynamics of non-uniform black strings. 

To calculate asymptotic charges such as mass and tension, one has to know the asymptotic behavior 
of the metrics. As we will see in the following sections, homogeneous perturbations (zero modes) appear 
at the second order, and they decay according to some power law in the asymptotic region r 3> r+. 
Suppressing all exponentially small corrections, which comes from Kaluza-Klein (KK) modes in the 
first- and second-orderperturbations, we find that the leading asymptotic behaviors of the metric 
functions for d > 5 are Qj 

a(r, z) ~ — 
V r 

/r,\ d—3 

\ r 

c(r,z) -d^ + Coot^) , (13) 
r V r J 



1 For d = 4, the asymptotic form of c(r, z) is exceptional: c(r, z) ~ C OQ (r + /r) ln(r/r+) 
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where Aoo, B^, Coo and C\ are constants. Prom the Einstein equations in the asymptotic region, 
these constants are related as 



Aoo + 2£oo + (d - 4)C 



0. 



(14) 



With these asymptotics, the mass per unit length M/L and relative tension n (and tension T) in the 
z-direction (z € [0, L]) for d > 5 are computed [32j[33] as 



M 

T 

n = 



a, 



„d-3 



2(2d - 5)Boo + 2(d -2)(d- 4)Coo + d - 
2{d - 4)5^ + 2{d - 4)^ + (d - 2)q d -z 



167rG d+ i 
LT _ 1 - 

~M ~ {d^ 2)[1 + 2(d - 4)Coo] + 2(2d - 5)B M + q 



2 + 



d-3 



d-3 ' 



(15) 
(16) 



where fi d _ 2 = 2tt {<1 ~ 1 ^ 12 /T[(d - l)/2] is the surface area of a unit (d — 2)-sphere. 

Entropy 5, temperature T, magnetic charge Q, and chemical potential <&h are given by 



m d-3 

T = ~A \ 

47rr + 



r=r+ 
~/„&+(<f-2)c\ 



/- e a ~ b 
//-<« 



i 4G d+1 + v 
Q _ (d- m d _ 2 r d + ' 3 
L 



r=r_|_ 



(d-3)/2 



16vrG d+ i 

<D H = ( d - 3)r d - 3 q^/ 2 [°° -^-{e a+2b -^ c ) 



(17) 
(18) 

(19) 
(20) 



r+ 



where ( • ) denotes an averaging in the z-direction. Here, some remarks on these quantities may be 
needed. Although the z-independence of the temperature (corresponding to the zeroth law) is not 
obvious from the expression (|17|) . it will be verified from the boundary conditions which state that a — b 
for all modes vanishes on the horizon except for the homogeneous mode. The magnetic charge (|19p 
is normalized so that we have Q/M — ► 1 in the extremal limit q — > 1. The chemical potential (|20l) 
can be calculated in two ways. One is to use the thermodynamical relation = (dM/dQ)s- F° r 
uniform black strings, this method is demonstrated in Appendix O The second way is based on the 
electric potential, by taking the dual field of the field 3~. We integrate out the dual field along the 
z-direction, by the KK reduction, and introduce the electric potential on the horizon by &h = -A.\ r=r+ , 
where *£F = dA. 

According to the correlated stability argument [231 124] . the dynamical stability of the background 
spacetime is closely related to its thermodynamical stability. In our case, the thermodynamical sta- 
bility is determined by the sign of a specific heat, given by (Appendix ICj) 



dM 
~dT 



m 



d-2T 



d-2 



,d-3 



[(d-2) 



,d-3i 



4G 



d+l 



-l + (d-2)q a 



(21) 



The specific heat Eq. pip is negative for a small charge, but it becomes positive at a critical charge 
Qgm, which is given parametrically by q = qcM, 



qGM 



(d - 2)V(*-3) 



(22) 
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Thus the system will be stable for Q > Qgm- As we will see below, the GL critical mode indeed 
disappears at this point, and the correlated stability is realized. 

Note that the choice of thermodynamic ensemble depends on crucially on whether we treat the 
charge as a parameter that has been fixed or as a parameter that can vary freely |34} [35l [TUl [36] . The 
above condition on the thermodynamic stability, Cq > 0, is for a canonical ensemble, in which the 
temperature and charge are kept fixed. If the solution is smeared, in a magnetic case this means the 
string is charged along a transverse direction, the charge itself can redistribute freely in the smeared 
direction, and we have to consider the stability in a grandcanonical ensemble in which the stability 
condition is specified by the specific heat and isothermal permittivity. 



3 Static linear perturbation : Xi$ 

3.1 Search for static mode 

We can specify the critical wavenumber k$ for a given value of q by solving the coupled ODEs fj 1 01 ) 
with suitable boundary conditions. Solving the EOMs in the asymptotic region (y ^> 1), we see that 
the perturbations resulting in regular solutions must have the following asymptotic behaviors: 

ai,o, ci, ~ e~ k °y. (23) 

From these asymptotics, it is understood that ai,o(y) and c\fi(y) are the KK modes, which are localized 
near the horizon. 

The event horizon is a regular singular point of the differential equations. Demanding the regularity 
of the perturbations on the horizon (y = 1), the following boundary conditions are required (for d = 5): 

, _ 2(kl - 6 + 16<? 2 - 18g 4 )ai - 3(fcg - 4 + 2Aq 2 - 20g 4 )ci 

a i,ol?/=i " 6(1 - g 2 ) 2 

-4 + kl + 16<? 2 - 12g 4 
2(1 - q 2 ) 



Ci,oly=i = 2 «i,o + —-t Cl fi- ( 24 ) 



Thus, these boundary conditions constrain the values of perturbations and their first derivatives. Since 
we are working on linear perturbations and the amplitude of perturbations is specified by e, we are 
free to fix the amplitude of a function at any point. Therefore, we set ci.o(l) = 1. The regularity 
at the horizon requires also &i,o(l) = a i,o(l)> corresponding to the zeroth law of thermodynamics, as 
mentioned in Sec. 12.31 

Now, we are ready to solve (I10p for each q with the boundary conditions (|23H and (I24p . For the 
numerical method, we integrate (|10p from the horizon by taking 0^0(1) and ko as shooting parameters. 
The shooting parameters resulting in regular solutions for the neutral case (q = 0) are given in Tables [2] 
for several dimensions (D = 6, 10 and 14). The charge dependence of the critical wavenumber is 
depicted in Fig. [Q We see that k$ decreases monotonically as the background charge increases, and 
then vanishes at the critical charge (|22p . above which the specific heat of the background solution is 
positive. This is a clear realization of the correlated stability. 



3.2 Transition 

The emergence of the static mode and its associated non-uniform states at Q < Qgm looks like a 
phase transition in condensed matter physics. H A simple example of phase transition is as follows. At 

2 We would like to thank Barak Kol, who pointed out this to us. 
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0.2 0.4 0.6 0.8 1 
Q/Qgm 



Figure 1: Charge dependence of GL critical wavenumber ko, normalized by the inverse of the hori- 
zon radius, for D = 5 ~ 15. We see that the critical wavenumbers decrease monotonically as the 
background charge increases. As expected from the correlated stability conjecture, these decreasing 
wavenumbers vanish at the critical charge Qgm, above which the background uniform black strings 
are thermodynamically stable. 

high temperatures, there is no order, and the "order" parameter ((f)) is zero. At a critical temperature, 
T c , order sets in so that, for temperatures below T c , (0) is nonzero. If ((f)) rises continuously from 
zero, the transition is second order and the specific heat has a jump discontinuity at the transition. 
[37] . The same thing holds for the present case, by just replacing (eft) and T by ko and Q, respectively, 
although the "order" (translational symmetry) sets in at Q > Qgm- The specific heat ([2T]) also has 
a discontinuity at Q = Qgm- This kind of phase transition can be well described by the Landau 
phenomenological mean-field theory. The mean-field theory, in which the free energy takes the form 
of / = r (4>) 2 + u (4>) A with u > 0, predicts a second-order phase transition with 

(4>)~{T C -Tf> (25) 

where /3 = 1/2 is called the critical exponent. It is interesting to see if the emergence of the static 
mode follows the line of this argument. The critical exponents extracted from the results presented in 
Fig. CD are 

/? = 0.55 (-0 = 6), = 0.51 (D = 7), (3 = 0.50 (D > 8), (26) 

These values are fairly close to the prediction of the Landau theory. (Note that we can also analyze 
the critical exponent for the mass per unit length, i.e. koM/2ir, as an another physically motivated 
quantity. Such analysis gives the same result as above because the contribution of q in M is negligibly 
small.) P in D < 7 deviates slightly from the prediction. The determination of fco near the GM point 
requires high numerical accuracy (in particular for smaller dimensions), and hence the derived critical 
exponent will contain numerical errors. Taking this into account, we are well justified in expecting 
that the critical exponent for kor + (or koM/2ir) is exactly f3 = 1/2 and that it will not depend on the 
dimensions. Moreover, it is an interesting issue to study the universality of this phenomenon. (See 
[38] for an earlier study of this kind of observation in other systems.) However, it is beyond the scope 
of this paper, and we leave it for future study. Instead, we extend the perturbative analysis to higher 
orders and identify the physical properties of charged non-uniform states. 
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4 Static perturbations at higher orders 



In this section, we discuss the structure of the EOMs, boundary conditions and asymptotics for higher- 
order perturbations. For simplicity, we set the dimension to d = 5 in the rest of this section, although 
the analyses are done for several spacetime dimensions. The dimensional dependence appears mainly 
in the asymptotic behavior for the zero modes at the second order. When we solve the ODEs at each 
order, a non-trivial task is to learn how to fix the integration constants. We will also discuss this point 
to some extent, in addition to other technical aspects. Readers who are interested only in the physical 
results can jump to Sec. Gf| 

4.1 Second-order perturbation 

At the second order 0(e 2 ), we have two independent modes, X 2 ^ and Xqj_. The former is a second- 
order counterpart of X\ o> which is straightforward to integrate. On the other hand, the latter is the 
lowest-order zero mode. This mode has a different structure of EOMs and thereby the asymptotic 
behaviors are different from those for KK modes. 

4.1.1 KK mode: X 2j0 

The EOMs for the X2 mode take the form of 



L [ 2l o a2,o + P [l] C2,o = S [ ^(X 1)0 ;k ), 

41 c 2,o + pl2] a2,o = S [ 2 2] (X lj0 ; ko), (27) 



similar to X\ t Q. Here, S 2 Q represents a source term, which is quadratic in Xio(y) and its first 
derivative. The function 62,0(2/) is written in terms of 122,0 (y), c 2,o(y) an d their first derivatives. 
The asymptotic behaviors in the far region are 

02,0, c 2 ,o ~ e~ 2k °y. (28) 
The horizon boundary conditions analogous to Eq. (|24p can be immediately obtained from ()2T|) : 

4,0(1) = 12 (1 1 g 2)2 [ 8 (~ 3 + 2k + 8 ^ " 5 9 4 ) a2 >° " 24 (" 1 + k + ~ 5( ? 4 ) C 2,0 

+2(-6 + k 2 + IQq 2 - 10g 4 )a 2 j0 + 3(8 + k 2 - 48<? 2 + 40g 4 )ai, ci, - 3(4 + 3k 2 - 64g 2 + 60g 4 )c 2 j( 

c 2 ,o(l) = 4(1 ^ 2 )2 [ 8 (" 1 + k o + V " 3 ^ 4 ) c 2,o + 8(-l + g 2 ) 2 a 2 , 

+4(-l + q 2 ) 2 a\ fl + (-8 + k 2 Q + Z2q 2 - 24g 4 )ai, ci, + (4 + 3A: 2 - 40g 2 + 36g 4 )c 2 ] , (29) 

where all functions are evaluated at y = 1. The regularity also requires 62,0(1) = <~2,o(l)- 

By taking ci2,o(l) and C2,o(l) as shooting parameters, one can easily integrate the coupled ODEs 
((27D with the boundary conditions ^2SJ) and 



3 A Mathematica notebook including numerics in this paper will be available at the author's website, 
http : //www. gravity. phys.waseda.ac.jp/~umpei. 
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4.1.2 Zero mode: Xq i 



The structure of EOMs for the Xq i mode is different from those for inhomogeneous modes. The 
starting point is to decouple 60,1 and co,i from 00,1 ■ @ The resultant equations for &o,i(y) an d co,i(y) 
take the form of 



1 



y 6 [5q*-3q 2 + 9y* + q 2 (l-12y 2 )]b' 



y 2 (-l + y 2 )(2g 4 -5gV + 3y 4 

-3g 2 y((? 2 - y 2 )(-2 + y 2 ) C ' li0 - 2q 2 [2q 2 + 3y 2 (-2 + y 2 )]6 0jl - 6g 2 (2g 2 - 6y 2 + 4y 4 )c 0il 

= ^cu (^1,0; &o) j 



c o,i 



y 2 (-l + y 2 )(2g 4 -5gV + 3y 4 ) 



y [3y 4 (-3 + Ay 2 ) + g 4 (-4 + by 2 ) + g 2 (lly 2 - 15y 4 )] c' 0>l 

-V(q 2 ~ 2/ 2 ) 2 (-2 + 3y 2 )b' 01 + [-4q 4 - 6y 6 + 2<rV (4 + y 2 )]6 0il + 2[6g 4 + 3y 6 - gV (12 + y 2 )]c ,i 

= S 2 \ (X lfi \ feo) , (30) 

where Sq^ represents a source term, which is quadratic in X±fi and its first derivative. 

Since the source terms decay exponentially in the asymptotic region, we can neglect them there. 
By the asymptotic expansion, we find 60,1 (y) an d CQ t i(y) have power-law asymptotics, which are given 
in Eq. fjl3[> . Here, noting that the constants in Eq. (| 13j) are 0(e 2 ), we write them as ~ e 2 b oc , 
Coo - e 2 Coo, C\ ~ e 2 ci (and A oo ~ e cioo)- From the regularity at the horizon, the horizon boundary 
conditions are 



b' 0il \ y =i = g 2 (6c ,i - a 2 - 2& ,i) + 



;2 

«, 



4 V 1 - q 2 

2\„_ 1 /1 „2w„2 1 r>;,2 



2\ 3 



ft m„2 



cb,ily=i = "2(1 - 3g 2 )c ,i + (1 - q 2 ){a( „ + 26f 



The equation for ao,i takes the form of 

1 



4 VI 
+ kl - 32q 2 + 24g 4 



12g^ cf o, (31) 



-«i,oci,o 



4(1 - g 2 ) 

4 - 3/c 2 - 40g 2 + 36g 4 2 



+ 



4(1 - g 2 ) 



to- (32) 



y 2 (y 2 - g 2 )(y 2 - 1) 



y[g^(2 + y 2 ) - 32/ 4 ]a' 0il + 3y(g 2 - y 2 )^ + 4g 2 (6 ,i - 3c ,i) = S\\(X lfi ; k ). (33) 



The asymptotic behavior is given by Eq. (|13|) again. The regularity at the horizon requires 
, , 1 



2(6 + k 2 - 16q 2 + 10g 4 )a 2 - 8(1 - q 2 )[(5q 2 - 3)6 ,1 



*0,lllf=l - 12 (l_g3)2 

+3(1 - 5g 2 )c ,i + 3(/cg - 8 + 48g 2 - 40g 4 )a lj0 ci i0 + 3(4 - 3k 2 - 64g 2 + 60g 4 )c? ]1 , (34) 



4 This procedure is peculiar to this charged system. For the neutral case (q — 0), it is possible to decouple 60,1 from 
ao.i and co.i, and then we can solve 60,1 ■ Since such a equation for feo,i contains only derivatives of 60,1, there is a shift 
symmetry and thereby we can integrate it without shooting |7|. However, because of new terms coming from gauge 
fields, it is hard to obtain a single equation for 60,1. 



9 



Since ao,i is decoupled from 60,1 an d co,i, all the values in the right-hand side (RHS) will be known 
after the integration of 60.1 and co.i- Furthermore, we notice that the above equation (133j) and the 
boundary condition (|34p contain only the first and second derivatives of ao,i(y). This means that 
ao,i has a shift symmetry, ao,i — ► ao,i + constant, corresponding to the gauge degree of freedom of 
redefining a time coordinate. Therefore, after integrating (|33|) with an arbitrary value of 00,1(1), one 
can use the shift symmetry in order to realize ao,i — > at the asymptotic region. No shooting is 
necessary. 

Finally, we discuss the integration of Eq. (|30p for 60,1 and cq,i- It seems that we have two constants 
to be fixed by asymptotic flatness, the horizon values of 60,1 an d oqx- However, we face the known 
puzzle (for the uncharged case) [7] 1181 ITT] : any choice of co,i(l) results in a regular solution satisfying 
asymptotic flatness. The physical quantities depend on the choice of co,i(l), and thereby the physical 
quantities cannot be fixed at this order. This apparent contradiction can be understood by the fact 
that the general conformal ansatz ([3]) has the freedom to "superpose" arbitrary 0(e 2 ) change in the 
mass on the top of the change that the non- uniformity induces I18j . The same argument holds 
for the present charged case: the equation (|30j) with the source term set to zero, i.e., homogeneous 
equation, has a solution that corresponds to the s-wave perturbation of a ci-dimensional charged black 
hole. The increase of mass and charge of the black hole is expressed as a solution to the homogeneous 
equations with the asymptotics (|13|) . Any multiple of such a solution to the homogeneous equation 
can be added to the solutions to the inhomogeneous equations. This is the origin of the ambiguity. 

Indeed, the horizon values of homogeneous solutions can be related to the change of physical 
quantities as follows. Consider a uniform black string given by setting X(r,z) to X{r) (X = a,b,c) 
in Eq. (|7]). Introducing a conformal radial coordinate 



p(r;r±) 



dr' 



Vf+(r';r + )f±(r';r- 



(35) 



the metric can be written as 



ds 2 



d+l 



-e 2a f + dt 2 + e 26 /- (dp 2 + dz 2 ) + e 2c r 2 dn 2 „ 2 . 



(36) 



In Eq. (|36p . r should be regarded as a function of p, which we denote by r = r(p, r±). The integration 
constant of Eq. (|35p can be fixed so that p = <-> r = r+. Because of the three functional degrees 
of freedom, there is a solution that changes mass and charge of the background black string. Such 
solutions with new parameters r± are given by the following homogeneous shifts of metric functions: 



X(r(p;r±)) -> X(r(p; r±)) + A x (p; r±,r'±), 



A a 



I In 
2 



f + (r(p;r'±);r' + ) 
f+(r(p;r±);r+) 



A; 



iln 
2 



f-(r(p; r'±);r'_) 
f-(r(p; r±);r_) 



A, 



In 



r(p\ r'±) 
r(p;r±) 



(37) 



where we write f±(r) as f±(r;r±) in order to show explicitly its dependence on the parameter r±. 
Corresponding to these shifts, we must also change the charge that appears in the RHS of the Einstein 
equation through the field strength. With all these transformation, the resultant solution is the charged 
uniform black strings with new parameter r^.. It may be useful to derive the relation between the 
horizon values of the shifts and parameters r^_. Taking the limit p — > in (|37p . we have 



lim A 



In I r -± 

r + 



+ In 



1-q 



,d-3 



lim A5 



iln 
2 



1 



Jd-3 



-.d-3 



lim A r 



In 



- ± ), (38) 
r+J 
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D 


0-1,0 


k 


«2,0 


C2,0 


«0,1 


&o,l 


Ooo 


boo 


Coo 


0-1,1 


h 


6 


-0.78 


1.27 


0.64 


-0.87 


0.63 


0.91 


-0.93 


0.65 


-0.37 


-0.38 


1.02 


10 


-1.54 


2.30 


2.48 


-1.53 


0.66 


1.46 


-4.06 


1.67 


0.14 


-0.44 


0.12 


14 


-2.15 


3.03 


5.16 


-2.14 


0.43 


1.90 


-10.0 


2.69 


0.51 


1.77 


-6.13 



Table 2: Shooting parameters (horizon values of metric perturbations and critical wavenumbers) 
and asymptotic quantities for the regular normalizable solutions with Q = 0. See section [J] for their 
definitions. We can use these results as initial-guess values for charged cases, Q 7^ 0. 

where q' = r'_/r', and we have used linip^o r(p; r±) = r+. 

At the practical level of calculations (perturbative expansions, etc.), it is easy to change 00,1(1), 
but is cumbersome to perform the necessary shift of the parameter r± which comes in through the 
field strength. If we forget the required shift of the parameter r± in the field strength and change only 
co,i(l) in our approach, the new background charge is automatically fixed and we have no freedom to 
perform an arbitrary shift of &o,i(l). In fact, as Eq. (|38p predicts, our numerical solutions confirm that 

(a ,i - 2&o,i + c ,i) y =i (39) 

is invariant under the shift of co,i(l), in which 60,1 (1) is automatically determined as a result of the 
shooting procedure. In the present paper, we want to construct the charged non-uniform string for a 
fixed charge and compare its physical quantities at the fixed charge. For this reason and to perform a 
clear analysis, we set co.i(l) = and solve Eq. ([3"0]) by taking 6o,i(l) as the shooting parameter. The 
shooting parameter resulting in a regular solution and the asymptotic constants for the neutral case 
are given in Table [2j As a consistency check, we will show later that the final result does not depend 
on the choice of co,i(l) (see Fig. [3|). 

4.2 Third-order perturbation: X\ y \ 

In the present approach, it is necessary to know the correction to the GL critical mode, fei, which 
appears at third order 0(e 3 ). At this order, we have two independent modes, -X'i 1 and A^o- However, 
since the equations for the latter mode X^^ do not contain k\, it is sufficient for our purpose to solve 
A M . 

The structure of EOMs and boundary conditions for A^i are similar to those for the A^o and 
A2 5 o modes. By decoupling £>i,i(y) from ai,i(y) and c\ t i{y), we obtain EOMs for ai,i(y) and ci ; i(y), 
taking the form of 

Lgoi,i + pWc M = S [ ^{X lt0 ,X 2>0 ,X 0jl ;k ,k 1 ) ) 

Lf Q c lyl + P^a lyl = sf\{X 1>0 ,X 2i0 ,X oX ,k ,k 1 ). (40) 

The function &i,i(y) can be written algebraically in terms of ai,i(y) and ci 1(2/). By solving (|40p in 
the asymptotic region, we obtain asymptotic behaviors, 

My), ci,i(y) ~ e~ k °y. (4i) 

From the regularity at the horizon, a\ 1 (1) and c' x 1 (1) can be written in terms of 0,1,1(1), ci i i(l) and 
the values and first derivatives of X\q at the horizon. 
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Which parameters should be determined by shooting at this order? According to [7], the different 
choice of k\ corresponds to a different "scheme" (gauge). We adopt the "standard" scheme in which 
we set ci iP (l) = for p > 0, and permit k\ to varyo Setting ci 5 i(l) equal to zero, we determine ai,i(l) 
and k\ by requiring normalizable modes. Namely, we integrate Eq. (|40p with the horizon regularity 
condition and the asymptotic damping ([IT]) by taking ai,i(l) and k\ as the shooting parameters, as 
we do at the first order. 



5 Physical properties of inhomogeneous charged strings 
5.1 Changes of physical quantities 

With the formulae prepared in Sec. 12.31 and perturbative quantities calculated in Sees. 13.11 and [H 
the increase of thermodynamical quantities due to the non-uniform deformation can be calculated. 
Because the asymptotic size of the transverse circle is not fixed {k\ / 0), we introduce variables that 
are invariant under rigid rescalings of the entire solution. Such invariant quantities can be obtained 
by multiplying K by suitable powers. The relative increase of the variables is given by 

o> 5M . , n ,5K 2 ^. 4 , 

/i M K 

- = — + d - 2 — = Sl e 2 + O f 4 
s b K 

£ ,« +((i -3)f = ^ + 0( e «), (42) 

where the second-order coefficients (for d > 5) are given by 

2(2d- 5)6 00 + 2(d-2)(d-4)c 00 , 

m = d-2 + q*~* + ( *~ 3) fc' 

j, , i a o\ , [»i,o + (d- 2)ci, ] 2 . /, n ,h 

si = &o,i + (d - 2)c ,i + — ^— - i— ^- + (d - 2) — , 

4 fe 

». kl 
n = a ,i - 00,1 - i— , 

fc 

tf 1 = (d-3)^. (43) 

«0 

The potential and relative tension are invariant by themselves, and their relative increase is given by 
®H n 

$ m = (d - 3) ^ {"o. 1 + 26 o,i - (<i - 2 )co,i + J [«i,o + 26i j0 - (d - 2)ci, ] 2 | , 

2(d - l)(d - 3)[(1 + 2g d - 3 )ft 00 + (d - A)q d - 3 Coo ] 
Ul {d - 2 + q d -3)[l + {d - 2)q d ~ 3 ] ' 1 j 



J It is also possible to set fei = in addition to ci jP = (p > 0). However, this requires a technically difficult 
procedure. In this case, co,i(l) at the second order should be dealt with as a shooting parameter so that the third-order 
perturbations converge at all. It means that the second- and third-order perturbations have to be solved at the same 
time. 
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The charge dependence of these quantities are depicted in Fig. [2] for D = 6, 10, 14. For D = 6 and 
10, one can observe that the change of mass /ii, being positive initially at Q = 0, decreases as the 
background charge increases. Then, [i\ becomes negative, e.g., at Q ~ 0.5Qgm for D = 6. The 
increase of entropy and charge, s± and #i, behaves in a similar way. The temperature ri also behaves 
in a similar way but in the opposite sign; the non-uniform black string, being "cooler" than the critical 
solution initially at Q = 0, becomes "hotter" as Q increases. Some remarkable properties are that 
does not change the sign and n\ changes its sign near the GM point where s\, Q\, t\ also cross 
the zero axis. The positive value of m has an interesting physical meaning. The uniform black string 
is not necessarily the state that has the largest tension (for a fixed background charge). This is not 
impossible for uncharged black stings. 

In Fig. [2j we have shown the results up to Q/Qgm — 0.95. The physical quantities in Eq. (f43l) 
contain the term of k\/ko, which diverges in the limit of Q — > Qgm, and thereby numerical accuracy 
becomes relatively worse near the GM point. To confirm the consistency of analysis and the accuracy 
of numerical integration, we use the Smarr formula as an independent check. 

According to Ref. |39} 140]. we find the following result: 

(D - 3)£ -{D- 2)TS - T + t {mat ^ = 0, (45) 

where £ and 8 are the mass and entropy density along the z-direction, respectively, and T is the 
tension. g( mai ) i s the mass density for matter fields evaluated on the horizon. It is given by the 
surface integral on the horizon at z, 

g{mat){z) = Jf^J/ S MN Pl ...P p ? MNPl --- Pp <PH(r + ,z), (46) 

where (j>u(r + ,z) is the density of electric potential at z. Taking the dual of Eq. and substituting 
it into the above formula, the mass density for the matter field is found to be £( ma *) = (Q/L)</>#. 
Integrating (|45p over the z direction yields the Smarr formula which holds even for the non-uniform 
black strings: 

{D - 3)M - {D - 2)TS -LT + Q<$> H = 0. (47) 

This formula can be also rewritten in terms of relative tension n as 

TS ( Q\ d-2-n 

-m + —i* >h [m)=^i^- (48) 

Now we expand all physical quantities as power series of parameter e: 

oo 

y = v 2p > ( Y = M > s > etc -)> ( 49 ) 

Substituting the expansion into Eq. (|48p . we obtain the following relation at 0(e 2 ), 

5S, + 8T X - 6M 1 + ^-g— 8m + {d , 7 3) ^f_ g0 ( ^-^—Sm - SM X + 6 Ql + 5* m ) = 0, (50) 
a — 2 — no (d—l)boTo \d — 2 — uq ) 

where we abbreviate bY\ = Yi/Yq. We can replace the quantity 5Y\ by the scale invariant quantity, 
e.g., 5 Mi — > fix, in Eq. (|50p . Therefore, Eq. (|50p is a constraint equation that the physical quantities 
must obey. Substituting our solutions into the Smarr formula, we have confirmed that the Smarr 
formula holds within a relative error of only a few percent. 
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Figure 2: The differences of physical quantities (|42ft and (|43]) between the non- uniform and the critical 
uniform black strings for each value of background charge Q for D = 6, 10, 14. For D = 6 and 10 
(< 13) one can see that all quantities except &hi change their signs. To obtain the results, we have 
adopted cq,i(1) = 0. 



5.2 Thermodynamical stability in microcanonical and canonical ensembles 

We are ready to study the stability of the inhomogeneous charged black strings. First, let us focus 
on the entropy comparison in a microcanonical ensemble, namely, the entropy difference between the 
non-uniform black string and a corresponding uniform black string of the same mass and charge. Such 
a difference is given by 



■S*nu - S\j 
<Snu 



(71 = Si 



= aie 2 + a 2 e 4 + 0(e 6 ), 
d - 2 + <5> 2 m (d 



0"2 



(d 
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01, 
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2 
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(51) 
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Q/Qgm 

Figure 3: The difference of entropy a% (free energy p<i) between the non-uniform and uniform black 
strings of the same mass (temperature) and charge for D = 6, 10, 14. Qi. cr and Qn jCr with the vertical 
lines show the locations of the critical charges. In each figure, we plot three results obtained by 
setting Co, i(l) = 0, +1, —1, although they cannot be distinguished from each other. For D < 13, the 
non-uniform string with Qi )Cr < Q < Qn )CT has larger entropy than the uniform string. For D > 14, 
the free energy difference is never positive, and the first critical charge does not exist. 



These equations can be derived by starting from the frame with 5K = and re-expressing the results 
in terms of the invariant quantities (|43j) . (See Appendix [D] for the details of calculation.) The first 
law of thermodynamics requires a\ = 0. The vanishing of o\ indeed holds with permissible numerical 
errors (a few percent relative errors) except near the GM point, e.g., Q > 0.95(5gm- Thus, the 
difference of entropy appears at 0(e 4 ). Note that we have used a\ = to simplify the expression of 
02- If 02 ^ s positive, it means that the non- uniform state has larger entropy than the uniform state 
and the non-uniform state is more stable. If 02 is negative, the uniform state is more favorable. For 
the uncharged case, the latter case is realized for D < 13, and the former case is for D > 14. 

The dependence of o<i on the background charge Q is depicted in Fig. [3] for D = 6, 10, 14. As 
expected from the sign change of fx±, observed in the previous section, the sign of 02 changes depending 
on the background charge. For D < 13, 02, being negative initially at Q = 0, increases as the charge 
increases, and it becomes positive at some critical charge Q = Qi jCr . However, the increase is not 
monotonic, rather, 02 has a peak. Increasing the charge furthermore toward the GM point, a 2 falls 
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off and becomes negative at a second critical charge Q = Qn lC r- The numerical value of the second 
critical charge corresponds to the zero-crossing point of n\ (or the second zero crossing of s±, fi\, etc.). 
This behavior near the GM point is general. For instance, o"2 for D = 10 shows similar behavior. 
A quantitative difference is that the first critical charges are achieved by smaller charges, while the 
second critical charge does not change so much. For D > 14, since a"2 is initially positive at Q = 0, 
the first critical charge that we observed for D < 13 is absent, and only the "second" critical charge 
exists. Above the critical charge, the non-uniform strings have less entropy, as we see for D < 13 near 
the GM point. 

The sign change of ai is an interesting phenomenon. In the present system, the charge works 
as a parameter that controls the order of transition between the uniform and non-uniform states. 
For Qi >cr < Q < Qn. cr in D < 13, the stable phase of the non-uniform state emerges, realizing the 
second-order transition from the uniform phase. In addition to this, a second remarkable feature is 
that the non- uniform states become less stable near the GM point {Q > Qii )Cr )- It is hard to provide 
a physical reason for these phenomena^. However, the second feature may be related to the fact that 
the non-uniform state itself vanishes at the GM point. 

Next, let us focus on the Helmholtz free energy [F = M — TS) in a canonical ensemble, namely, 
the comparison of free energy between the non-uniform black string and a uniform black string of 
same temperature and charge. Such a difference is given by 

FNU " FU =p 2 e± + 0{e% (52) 



(d-3)(l-*% ) f (g-ljjlrogi 
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Sl+ l-(d-2)^ Tl 



(53) 



where we have used the first law of thermodynamics (Appendix [D]). Again the difference appears at 

0(6 4 ). 

The results are shown in Fig. [3] for D = 6, 10, 14. For the uncharged case, p2 is positive (for 
D < 12), and thus the uniform state is more favorable due to the smaller free energy. As expected 
from the results of entropy comparison, the sign of p2 changes as the charge increases. However, in this 
canonical ensemble, the sign changes only once. p2 tends to decease monotonically, and the critical 
charge is almost exactly the same as Qi iCr for the entropy comparison. It is remarkable that the critical 
charge does not exist for D > 13 and p2 is always negative. This result can be understood with the 
general behavior of p2 discussed above and with the fact that p2 for Q = becomes negative for 
D > 13. Therefore, the non-uniform phase in the canonical ensemble for D > 13 is stable irrespective 
of the charge. 

Finally, it is worthwhile to note that in Fig. [3] we have also shown the results obtained by taking 
co,i(l) = +1,-1, in addition to our standard choice co,i(l) = 0. The three lines overlap completely, 
and we cannot distinguish between them. As we have discussed in Sec. 14.1.21 there is a subtlety in 
the choice of co,i(l), but the present result means that our analysis is consistent and co 5 i(l) does not 
affect the final results, as expected. 



The global thermodynamic consideration, in which the entropy is compared between a uniform black string and 
a localized black hole of same mass (and charge), provides us an intuitive guess for the critical dimension/charge (e.g. 
8, 29]). The present system, however, is out of the global argument since the corresponding localized black hole solution, 
if any, cannot couple to the gauge field of the same type. 
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5.3 Comparison with the criticality of uncharged system 




5 7 9 11 13 15 5 7 9 11 13 15 

D D 

Figure 4: Dimensional dependence of physical quantities, e.g., si,//i, etc., in addition to the entropy- 
difference o"2 and free energy difference p2, for uncharged black strings (Q = 0). s± and jjL\ become 
negative at D = 14 and t\ becomes positive at D = 13. n\ does not change its sign, oi and p2 
change their sign at D « 13.5 and D « 12.5, respectively [HI [10] . These results suggest that the phase 
transition from the uniform to non-uniform black strings in a microcanonical (canonical) ensemble is 
first order for 5 < D < 13 (5 < D < 12) and higher order for D > 14 (D > 13). 

As we have already mentioned, we know that the stability of uncharged non-uniform string depends 
on the spacetime dimensions. (The number of extended transverse dimensions on tori T p is not 
irrelevant [11], and so we consider co-dimension one case of black string.) It is of great interest to 
compare the transition caused by the charge with the dimension-dependent transition of uncharged 
strings. In Fig.HJ we show the dimensional dependence of physical quantities of uncharged non- uniform 
string. We also plot the entropy difference and free-energy difference in Fig. [U 

Compared with Figs. [2]and[3l we see large differences in their qualitative behavior. The common 
behavior is that s±, [i\,T\ change their sign. For the charged case, the zero crossing of s±, fix, t\ occurs 
almost simultaneously, in contrast to the uncharged case. Because of this feature, the critical charges 
in the canonical ensemble and the microcanonical ensemble are almost identical, while the critical 
dimensions for the two ensembles are different. A second striking disparity is that the charged system 
has a second critical charge. This would be related to the fact that the GL static mode vanishes in 
the limit of Q — > Qgm while the static mode survives in the limit of D > 1 [9] . 

6 Summary and discussion 

We have constructed the perturbatively non-uniform black strings (NUBS) coupled with the gauge 
field. At the first order of perturbation, we have confirmed the realization of the correlated stability, 
i.e., the vanishing of the Gregory-Laflamme mode at the point where the background uniform black 
string (UBS) becomes thermodynamically stable. We have pointed out that the emergence/vanishing 
of the static mode resembles phase transitions. In fact, we have identified the critical exponent for 
the wavenumber of the static mode and found that the critical exponent is very close to (3 = 1/2, 
which means a second-order phase transition. However, this result does not necessarily mean that 
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the charged non-uniform black strings will arise from the uniform ones as a result of a second-order 
transition. This is because the black hole system is different in several aspects from the ordinary 
system in condensed matter (7j, and, in addition, ko does not directly specify the state after the 
transition, but specifies only the point where such a state emanates. 

The analysis of higher-order static perturbations allows us to study the charged non-uniform state 
with its (thermodynamic) stability. Our analysis shows that the charged non-uniform strings can have 
larger entropy than the corresponding uniform string in the appropriate range of charge (Fig. [3]). This 
means that by adjusting the background charge the smooth transition (evolution) from the uniform 
state to non-uniform state becomes possible. This is in contrast to the uncharged system: we have 
known that uncharged black strings can have larger entropy (and smaller free energy) than that of 
uniform ones of the same mass (temperature) only in large spacetime dimensions, D > 14 (D > 13). 

For the unchanged system, the phases of uniform and non-uniform states can be described by an 
analogous method of the Landau theory [T2] , where the external non-ordering field is the inverse 
of temperature. In the present case, the charge is a new non-ordering field, and because of this, phase 
diagrams should become multidimensional. Such a system can be also described by the Landau theory. 
Consider the following Landau free energy: 

/ ~ u 2 tf + u±(f> 4 + ■■■, (54) 

where the higher order terms maintaining stability are abbreviated. U2 is a function of the non- 
ordering field, for instance, u 2 oc (T™ 1 — TZ^) near the critical temperature in the uncharged case 
[111 112] . If «4 is positive, a second-order transition is predicted, while there is a first-order transition 
if U4 is negative. Like the normal-to-superfluid transition in liquid helium mixtures, U4 can depend on 
non-ordering variables and its sign can vary. For the two non-ordering fields, the so-called tricritical 
point appears, at which three phases of disordered state and ordered states with first- and second- 
order transitions meet. Clearly, our charged system fits this argument. U4 depends on charge (and 
temperature). It is an interesting issue to extend and refine the argument, based on [HI [12] • In 
the following, instead of pursuing its precise implementation, we discuss the (M, n) phase diagram 
gathering all the information that we know. 

A mass-tension (M, n) diagram is commonly used to classify the black objects in the KK compact- 
ification. We project out the degree of freedom of charge on the (M, ra)-plane and consider the phase 
on this plane. In Fig. [5l we illustrate a (M, n) diagram for a charged system (D < 13). Although 
L oc 1/ko diverges in the GM limit {Q = Qgm)> the mass per unit length remains finite, so that we 
adopt it in the figure. The sequence of the marginally stable charged uniform black strings (UBS) is 
given by Eqs. (115D and (|16j) . and the non-uniform branch can be read off from Fig. [2j Corresponding 
to Hi < in <5i 5 cr < Q < Qn,cr> the branch of NUBS declines in the region. However, the corrections 
to the mass and tension become > 0, n\ > 0, and then the branch of NUBS ascends near the GM 
point. These behaviors show that the phase space of black hole and string has a very rich structure. 

Basing our ideas on this phase structure, let us speculate about a possible dynamical evolution 
of sufficiently charged black strings with a 2 > 0. First of all, we adopt the plausible assumption 
that the conservation of magnetic charge prevents pinching-off uniform black strings/branes. (This 
is also expected from the fact that the black hole cannot couple to the same type of form field.) 
Thus a possible scenario is that once the uniform black string reaches the GL point, via radiation, 
etc., the unstable critical string will transit to a non- uniform black string adiabatically if the charge 
is appropriately chosen so that the non-uniform branch is less massive and entropically favored. At 
present, the dynamical stability of non-uniform black stings is not known. However, it is possible 
that a shrinking part of the wavy (apparent) horizon will locally evolve into the extremality, and the 
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Figure 5: An illustration of charged non-uniform black strings (NUBS) in a (M, n) plane, by projecting 
the charge axis on (M, n) plane (D < 13). The branch of uncharged NUBS and black hole are based 
on [20]. Since the GL wavelength, and hence the bare mass, diverges in the limit Q — > Qgm, we 
plot the mass per unit length, M/L, where L is the critical GL wavelength L = 2n/K. The vertical 
dot-dashed lines represent uniform black strings with fixed changes, and the filled circles shows the 
sequence of critical GL points of a charged system. Since there is no static mode for Q > Qgm, the 
sequence of the critical GL point terminates at Q = Qgm- The lines with arrows show the direction 
in which the branch of charged NUBS emanates. The direction of arrows can be read off from Fig. [2j 



stabilization works locally, making the whole black string stabilized in the end. Therefore, a wavy 
black string would be a possible end state in such a dynamical evolution. To identify the sequence of 
evolution, it will be useful to construct charged non-uniform black strings non-perturbatively. Since 
there is no corresponding branch of black holes carrying the same type of charge, we have no reason 
to expect that the branch of black string remains near the uncharged one even for a weak charge. 
Therefore, it is interesting to know how the branch deviates from that of an uncharged string at large 
non- uniformity. 

So far, we have discussed the phase space in Q < Qgm where the static modes exist. However, the 
phase space in Q > Qgm may be also interesting to explore. An argument proposed by Horowitz and 
Maeda [H] is that (stable) non-uniform black strings would exist near the extremality. Such a solution 
has not been discovered so far, and it will be interesting to know how the branches of near-extremal 
non- uniform black strings emerge and how they are embedded in the diagram. However, we should 
recall here that the scenario of such new branch is possible only when cycles on the horizon cannot 
shrink to zero size. So if it can pinch off, a near-extremal non-uniform black string is not an inevitable 
result. 

Finally, let us comment on possible extensions of the work. For simplification, we have not taken 
into account the dilaton, but it will be straightforward to include the dilaton and other types of 
gauge fields. As discussed above, it is important to figure out the phase structure in fully non- 
linear regimes in the present system. This problem will be addressed in our forthcoming paper. The 
dynamical evolution of GL instability in a charged system will be the most challenging problem. To 
see what happens dynamically when the charge crosses the critical value will be quite fascinating. The 
dynamics of a string's evolution may change drastically at the critical charge. This can be achieved 
by a straightforward extension of the analysis in [1], and the result will also provide an expectation 
for the dynamical evolution of uncharged black strings above the critical dimension {D > 14). 
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A Einstein equations in static axisymmetric spacetimes 



We present the Ricci tensor components for the general metric ansatz ([3]). Here, we denote the partial 
derivatives by d r A = A' and d z A = A. In addition, we specify the coordinates of the angular part by 
dVt 2 l _ 2 = d6f + sin 2 Q\dQ\ + • • • + sin 2 Q\ ■ ■ • sin 2 9 c i-sd8 2 ,_ 2 - The Ricci tensor is calculated as follows: 



Rtt : 



2A-2B-2H I 2H 



e zn A 



{d-2)H> c2H& 



A + {d- 2)C 



+ A' 



A' - H' + (d - 2) ( C + 



+ e 2H A + A'' 



A + (d-2)C + A' (H' - A' + B') 



+(d - 2)C [H f - C" 



+ B' 
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H' — A' — (d- 2) [C' + - 

_ A[_ _ (d-3) _ ^ 
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R rz = (B- A) A' + A + (d-2)C B' + (d- 2) ( B - C ) ( C + 



1 



A' - (d - 2)C', 



(55) 



and the other angular components are given by R^ = R e e 2 2 = ■ ■ ■ = 2 . For the magnetic field given 
by F = Qsd-2, the non-zero components of the RHS of the Einstein equation ([2]), which we denote 
by Sfw, are 



Q^ c*r qz 

D t — D r — D z 



d-3 



a 1 



Si d — 1 r 2(d-2) e 2(d-2)C 

B Other form field components 



2(d-1) r 2(d-2) e 2(d-2)C" 

o 2 

'i = l,2,...,d-2). 



(56) 



In this paper, we focus only on the angular parts of the form field F. In d = 4 and d = 5, however, 
the form field can have other components, corresponding to an electric field. In this section, we briefly 
discuss these components compatible to the conformal ansatz (|3|). 



20 



For d = 5, F trz component can exist in addition to (|3J). We put this part as 

F e = E(r, z) U3t A u r A u z , (57) 
where Ut = e A dt, u) r = e B+H dr and oj z = e B dz. Substituting this into EOM ([2|), we obtain a solution 

where Q e is an integration constant. The solution (|58p is the general solution of electric field, which 
we ignored in our analysis. Note that since the 3-form is self-dual in 6 dimensions, the electric and 
magnetic parts contribute to the Einstein equation in a similar manner: if we consider the electric and 
magnetic part simultaneously, the source term of the Einstein equation is obtained just by replacing 
the magnetic charge Q 2 in (|56p by Q 2 + Q 2 . In other words, the electric field in d = 5 can be analyzed 
in parallel with the magnetic one, and therefore the self-duality does not affect the final results. 
For d = 4, the electric components in addition to ([!]) are possible and they are given by 

F e = E r {r, z) LO t A uj r + E z (r, z) oj t /\uj z + E rz (r, z) oj t A u z . (59) 

From the EOMs, we find E rz = Q r z/( r2 e 2C+A ), where Q rz is an integration constant. This is the 
general solution of r-z component. The other components are given by solving the EOM and Bianchi 
identities: 

(d r + 2r~ 1 )(e B+2C E r ) + e H d z (e B+2C E z ) = 0, 

d r {e A+B E z )-e H d z {e A+B E r ) = 0. (60) 



C Specific heat and chemical potential 

We derive the formula for the specific heat used in Sec. 12.31 First of all, note that every thermody- 
namical quantity of uniform black strings (a = b = c = 0) can be written as functions of r + and q. 
From Eq. (|19p . we have 



1 

^2" 



167rGd + i 



2/(d-3) 



(d- i)m d _ 2 

Substituting (f6T]) into the mass formula (fT5j) or the product of S and T, one finds 



(61) 



r ,_ 3 = WnG d+1 M 1] / 4(d-2)Q 2 



m d+ i 2{d - 2) I V (d- l) 2 M 2 



Let us first discuss the specific heat. Utilizing (|6ip and (|62p . the RHS of the temperature (|17h can 
be written in terms of Q and ST, i.e., T = T(Q, ST). Differentiating this equation with respect to 5, 
we easily find (dS '/ '8T)q, and the specific heat is evaluated by 

dM\ ^(dS 

Q 
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which yields the result ([2TJ) . 

The method explained above is the most straightforward for obtaining thermodynamic quantities, 
e.g., the specific heat. However, it is sometimes difficult to get an explicit expression of a thermo- 
dynamic quantity in terms of relevant thermodynamic variables. In such cases, it is necessary to 
perform rather messy calculations. Instead of such a direct method, we explain a systematic method 
to calculate thermodynamic variables which have parametric representations. 

As an illustration, we begin with the calculation of chemical potential. The chemical potential is 
given by 

8M\ fdM\ fdr+\ ( 8M\ ( dq \ 



*» s {M)r^)M) s + y^)rM) s - (64) 

To calculate (dr + /dQ)s and (dq/dQ)s, the charge Q(r + ,q) is brought into the form of Q = Q(r + ,q(S, r+)) 
= Q(r + (S,q),q). Differentiating these expressions with respect to r + and q as 

dQ\ fdQ\ fdQ\ ( dq 



dr +/S \ dr +Jq \ d( lJr+\ dr 



+ 



S).-(£H£).(£>- 



and substituting them into (|64p . we finally obtain <&h — Q^ d 3 ^ 2 . By replacing variables suitably in 
the above calculation of chemical potential, we can easily calculate the specific heat Cq. 



D Entropy and free energy difference 



We use the expansion (|49h again for all relevant physical quantities. The entropy difference between 
the non-uniform and uniform black strings for same mass and charge is given by 



SnIJ - S\j 
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Si _ AffW 
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+ 



S 2 AS® Si ( AS^ \ 
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Dq DQ Dq Dq V Dq J 



e 4 + 0(e 6 ), (66) 



where AS = ^ P =i AS^ p 'e 2p is the difference of entropy between the critical uniform black string and 
the uniform black string that has the same mass and charge as the non-uniform black string. 

Si /So and S2/S0 in ([66]) are easily computed from the first law for a fixed asymptotic length of 
the circle, dM = TdS + &ndQ. Integrating the first law with the expansion, we obtain 



M x =T S 1 + $ H0 Qi, 

M 2 = T S 2 + -TiSi + $ m Q2 + \$mQi. 



(67) 



Furthermore, from formulae (|15|) . (j!7|) . (|18|) . and (|19|). we can find that the mass of uniform black 
string Mq can be written in two ways: 



M n 



d - 2 + q 
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(68) 
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The first equality can be also confirmed by the Smarr formula. From (|67p and (|68p . we can write 
Si/S and S 2 /S in terms of M x /M , M 2 /M , Qi/Q , Q2/Q0, and Ti/T . 

Here, we focus on the quantity AS/ So in (|66p . As mentioned before, AS 1 is the entropy change of 
the uniform black string due to the change of mass and charge. Therefore, AS can be computed by 



AS 



OS \ 
dM , 



AM + 



dS_ 
dQ 



AQ 



M 



where 5, M, and Q are for uniform strings (a = 6 = c = 0). One can obtain the explicit dependence 
of S on M and Q by substituting the expressions of q and r + , Eqs. (jBTj) and (|62|) . into Eq. (fT5l) . But, 
it seems to be easiest to calculate ([69]) as follows. The partial derivatives of S with respect to M and 
Q in Eq. (f69|) are given by 
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(70) 



From Eq. ([19]), r + and q are rewritten as r + = r + (Q, q) and q = q(Q, r + ), so M(r + , g) can be written 
as M = M(r+(Q,q),q) = M(r+,q(Q,r + )). Differentiating these relations with respect to r + and q, 



we obtain the coefficient in Eq. IjTOjl . Then, by setting AM = Mie z + M 2 e 4 and AQ = Q x t l + Q2e 
obtain AS = AS^e 2 + AS^e 4 . 

We can find the leading term to be the same as the one for the non-uniform black string: 
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(71) 



where we again abbreviate 8Y\ = Yi/Yq (Y = S, M, Q). Hence the term of order 0{e 2 ) vanishes 
in the entropy difference and the leading order correction comes from <3(e 4 ). Note that M 2 /M and 
Q2/Q0 cancel out in the final expression as expected from the cancellation at the order 0(e 2 ). In the 
end the entropy difference for the same mass and charge is given by 
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(72) 



where we have used the relation (|7ip to simplify the result. Last, by replacing 5Y\ for invariant 
quantities, we have Eq. (f5Tj) . 

The derivation of the formula for the free-energy comparison is possible in a similar way to that 
of entropy. First of all, the free energy [F = M — TS) of a uniform solution can be written in two 
ways from Eqs. CES}, (ft?]) . (fT8|) . and (fT9l): 
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The difference of free energy between uniform and non-uniform black strings is given by 



F + AT 



1. 



(74) 



where AF is the free-energy change of a uniform black string due to the change of temperature and 
charge, which we denote by AT and AQ, respectively. From the definition of the free energy and the 
first law, we obtain dT = — SdT + <E>#dQ. Integrating this by using the expansion (fl9j) . 



F 1 = -S T 1 + $ H0 Q 1 , 

F 2 = -S T 2 - IsxTi + § m Q2 + \*hiQi- 



(75) 



Using (731) and (JTSJ), F X /F G and T 2 /T are written in terms of Ti/T , T 2 /T , Qi/Qo, Q2/Q0, ^m/^m, 
and Si/Sq. 

Now, we focus on the term of AF/Fq in (|74p . A necessary correction of the free energy to the 
uniform string is 
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where the differential operators acting on F are evaluated by using the following relation: 
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Here, for example, the coefficient (dr + /dQ) T is computed as follows. From Eq. (fTTj) . r + and q 
are rewritten as r + = r + (To,q) and q = g(To,r + ), and the charge Q(r + ,q) can be written as Q = 
Q(r + (To,q),q) = Q(r + ,q(To,r + )). Differentiating these relations with respect to r + and q, we obtain 
the coefficient. 

We find that the correction of 0(e 2 ) vanishes as for the entropy difference. 
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